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ABSTRACT 

We present estimates of the basic physical properties - including size and albedo - of the extreme Centaur 2013 AZeo. 
These properties have been derived from optical and thermal infrared measurements. Our optical measurements revealed 
a likely full period of « 9.4 h with a shallow amplitude of 4.5%. By combining optical brightness information and thermal 
emission data, we are able to derive a diameter of 62.3 ± 5.3 km and a geometric albedo of 2.9% - corresponding to an 
extremely dark surface. Additionally, our finding of > 50 Jm“^K“^s“^'^^ for the thermal inertia is also noticeably for 
objects in such a distance. The results of dynamical simulations yield an unstable orbit, with a 50% probability that 
the target will be ejected from the Solar System within 700,000 years. The current orbit of this object as well as its 
instability could imply a pristine cometary surface. This possibility is in agreement with the observed low geometric 
albedo and red photometric colour indices for the object, which are a good match for the surface of a dormant comet 
- as would be expected for a long-period cometary body approaching perihelion. Despite the fact it was approaching 
ever closer to the Sun, however, the object exhibited star-like profiles in each of our observations, lacking any sign of 
cometary activity. By the albedo, 2013 AZgo is a candidate for the darkest body among the known TNOs. 

Key words. Kuiper belt objects: 2013 AZeo - Radiation mechanisms: thermal - Techniques: photometric 


1. Introduction 

The object 2013 AZgo is a recently discovered extreme 
Centaur, moving on an eccentric orbit with e « 0.992 and a 
perihelion distance of g ~ 7.9 AU. As a result, 2013 AZgo is 
among the TNOs with the largest known aphelion distance 
at Ri 1950 AU. 2013 AZgo may be cl assified as a Centaur , 
based on its perihelion distance (|Horner et al.l . l2003li . 
However, due to its large semimajor axis, it could equally 
be co nsidered to be a scattered disk object (iGladman et ah . 
2008). Its Tisserand parameter (|Duncan. Levison fc Donea . 
2004) w.r .t. Jupiter is Tj = 3.47 which is typical for 
Centaurs (|Horner et al.l . l2004allbh and differs from that of 
Jupiter family comets (2 < Tj < 3) and especially for from 
that o f Damocolids and Halley-type comets (Tj < 2, see 
iJewittl . l2005l) that exhibit cometary dynamics. 

In order to recover the basic physical and surface char¬ 
acteristics of this object, we need measurements both in the 


visual and in the thermal infrared range. Optical data can 
yield information about the intrinsic colours, the absolute 
brightness, rotational period, shape and surface homogene¬ 
ity of the object, while thermal observations aid us to decide 
whether we see a “large but dim” or a “small but bright” 
surface. For thi s latter purpose, Herschel Space Observatory 
(|Pilbratt et al.l . l201flh is an ideal instrument since the ex¬ 
pected peak of the thermal emissio n is close to the shorte st 
wavelengths of its PACS detector (iPoglitsch et al.l . |2010D . 

In our current analysis, we follow the same method¬ 
ology as pr esented in our pr evious study of the Centaur 
2012DR 30 l|Kiss et al.1 . [201^, another object moving on 
a similar orbit to 2013 AZgg. The structure of this pa¬ 
per is as follows. In Sec. [U we describe our observa¬ 
tions, including the detection of thermal emission by 
Herschel/PACS, optical photometry by the IAC-80 tele¬ 
scope (Teide Observatory, Tenerife, Spain), optical re¬ 
flectance by the Gran Telescopic Canarias (GTC, Roque 
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de los Muchachos Observatory, La Palma, Spain) and near- 
infrared photometry by the William Herschel Telescope 
(WHT, Roque de los Muchachos Observatory, La Palma, 
Spain). In Sec. [31 we derive the basic physical properties 
of the object by applying well understood thermophysical 
models. The dynamics of 2013 AZeo are then discussed in 
Sec. m Finally, our results are summarized in Sec. |SJ 


Table 2. Thermal fluxes of 2013 AZeo derived from our 
Herschel measurements. 


Band 

A 

Flux 

B 

70 pm 

32.5 ± 2.2 mJy 

G 

100 pm 

23.0 ± 2.8 mJy 

R 

160 pm 

15.9 ±4.5 mJy 


2. Observations and data reduction 


2.1. Thermal observations and flux estimations 

Thermal infrared images have been acquired with the 
Photoconducto r Array Camera and Spectrometer (PACS 
iPoglitsch et al.L 2010ll camera o f the Herschel Space 
Observatory ( Pilbratt et al.l . l2010ll in two series, each of 
1.3 hour duration. As we summarize in Table. [H these 
two measurement cycles were separated by more than 
4 hours, allowing the target object 2013 AZeo to move, 
but still be in the same field of view. This type of 
data collection has almost exclusively been employed in 
the “TN Os are Cool!” Open Time Key Programme o f 
Herschel (|Miiller et al.l . l2009l 120101 : IVilenius et al.l . l2012ll . 
For both series of measurements, we used both the blue/red 
(70/160/rm) and green/red (100/160/rm) channel combi¬ 
nations. This scheme allowed us to use the second series of 
images as a background for the first series of images (and 
vice versa) in order to eliminate the systematic effects of 
the strong thermal background. This type of data acqui¬ 
sition and the respective reduction scheme were described 
in our former w orks related to both the “TNOs are C ool!” 
project (see e.g. IVilenius et all 120121: Pal et al.l. 201211 and 
subsequent measurements (see e.g. iKiss et al. L I2013D . 

Unfortunately, the astrometric uncertainties of 
2013 AZgo were relatively large at the time of Herschel 
observations, due to its rather recent discovery. Thus, 
the apparent position of the object was slightly (fsi 29") 
off from the image center, which also implied that the 
double-differential photometric method (iKiss et al.l . l2014ll 
yielded larger photometric uncertainties. In addition, 
shortly before the Herschel observations, on February 16, 
2013 (at OD-137^1), one half of the red (160 /rm) channel 
pixel array became faulty. Hence, only the images from the 
first visit were sufficient to obtain fluxes at 160 /rm and it 
was not possible to create double-differential maps in this 
channel. 

Raw Herschel/P ACS data h as been processed in the 
HIPE env ironment with custom scripts de¬ 

scribed in iKiss et al.l (|2014[) . The double-differential maps 
were creat ed and analyzed using the FITSH package 
(EIIIqiI. The resulting images are displayed in Fig. |T] 
Photometry on the individual as well as on the combined 
double-differential images were performed by using aper¬ 
ture photometry where the fluxes were corrected by the 
respective growth curve functions. Photometric uncertain¬ 
ties were estimated by involving artificial source implan¬ 
tation in a Mon t e-Carl o fashion (ISantos- Sanz et l2012t 
iMommert et al.l . l2012t iKiss et al.L 201^ This method 
works both for the double-differential images (blue and 
green channels) as well as on the individual maps (here, 
the red channel). 


^ http://herschel.esac.esa.int/Docs/Herschel/ 
/html/ch03s02.html#sec3:DeadMat 


70 /rm 100 pm 160 pm 



aR.4. (orcssc) AR.A. Corcsec) 


Fig. 1. Image stamps showing the Herschel/PACS maps of 
2013 AZgo in the 70 pm (blue), 100 pm (green), and 160 pm 
(red) channels. Each stamp covers an area of 64" x 64", 
while the tick marks on the axes show the relative positions 
in pixels. The effective beam size (i.e. the circle with a di¬ 
ameter corresponding to the full width at half maximum) is 
also displayed in the lower-left corners of the stamps. The 
upper panels show the stamps directly combined from the 
individual frames where the lower stamps are obtained us¬ 
ing the double-differential method. Due to the failure of the 
half of the red channel and the astrometric uncertainties, 
the second visit is nearly unusable in red, hence double¬ 
differential maps cannot be created. 


Based on the individual images, we obtained thermal 
fluxes of 36.6 ± 2.9, 25.2 ± 3.7 and 15.9 ± 4.5 mJy in the 
blue, green and red wavelengths, respectively. By involv¬ 
ing the double-differential maps, we derived 32.5 ± 2.2 and 
23.0 ± 2.8 mJy in the blue and green regimes. Due to the 
lower level of confusion noise (see Fig. [T]), the accuracy 
of the latter series of fluxes is better. Therefore, for fur¬ 
ther modelling we adopt the double-differential fluxes for 
blue and green. Thermal fluxes should undergo colour cor- 
rection according to the temperatures of the bodies (see 
IPoglitsch et al for the respective coefficients). Since 

the subsolar temperature of 2013 AZgo is around 110 K, the 
colour correction is negligible (less than a percent) in blue 
and green while it is -1-4% in red. Our reported fluxes con¬ 
sider the respective colour correction factors. In addition, 
in the error estimation of the aforementioned fluxes, we 
included the 5% systematic e rror f or the absolute flux cali¬ 
bration as well ([Balog et al.l . I20l3l . The summary of these 
thermal fluxes are reported in Table [3] 
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Table 1. Summary of Herschel observations of 2013 AZeo, obtained in the DDT program DDT_ckiss-3. The columns are: 
i) visit; ii) observation identifier; hi) date and time; iv) duration; v) filters configuration; vi) scan angle direction with 
respect to the detector array. 


Visit 

OBSID 

Date & Time 

Duration 

Filters 

Scan angle 



(UT) 

(s) 

(pm/pm) 

(deg) 


1342268974 

2013-03-31 18:10:51 

1132 

70/160 

70 

Visit-1 

1342268975 

2013-03-31 18:30:46 

1132 

70/160 

110 


1342268976 

2013-03-31 18:50:41 

1132 

100/160 

70 


1342268977 

2013-03-31 19:10:36 

1132 

100/160 

110 


1342268990 

2013-03-31 23:46:55 

1132 

70/160 

110 

Visit-2 

1342268991 

2013-04-01 00:06:50 

1132 

70/160 

70 


1342268992 

2013-04-01 00:26:45 

1132 

100/160 

110 


1342268993 

2013-04-01 00:46:40 

1132 

100/160 

70 



Fig. 2. Folded optical light curves of 2013 AZgo using 
photometric data taken on 6 subsequent nights of 2013 
November 4-9. Note that the folding frequency is related 
to the preferred double-peaked solution, n = (5.11/2) d“^. 
See text for further details. 


2.2. Optical photometry 

Since 2013 AZgo has recently been discovered, one of our 
goals was to obtain precise photometric time series for this 
object in order to estimate both the absolute magnitudes 
(in various passbands) and the rotational period based on 
light curve variability. For these purposes, we used the 
lAC-80 telescope located at Teide Observatory, Tenerife. 
During our observing runs, we used the CAMELOT cam¬ 
era, equipped with a CCD-E2V detector of 2k x 2k, with 
a pixel scale of « 0.3" providing a field of view of roughly 
10' X 10'. Time series were gathered for several hours on the 
nights of 2013 November 4-9 using Sloan g', r' and i' filter 
sets. Since 2013 AZgo can currently be found at the edge of 
the spring Sloan field, several dozen stars with accurate ref¬ 
erence magnitudes were available on each image. The night 
conditions were photometric on 5th, 7th and 9th of 2013 
November where the individual photometric uncertainties 



Wavelength (A) 

Fig. 3. Reflectance spectrum of 2013 AZgo, taken with the 
OSIRIS spectrometer on the GTC in January, 2014. This 
spectrum is normalized to unity at A = 5500 A. 


were nearly constant and varied between 0.04 and 0.05 mags 
in Sloan r' band. The conditions were worse on the three 
another nights (4th, 6th and 8th) when the formal uncer¬ 
tainties scattered in the range of 0.04 — 0.08 indicating the 
variable transparency of the sky (which was also notable 
during the observations). The worst conditions were on the 
first night where some of the measurements had a formal 
uncertainty of 0.09. 

The scientific images were analyzed with the standard 
calibration, source extraction, astrometry, cro ss matchin g 
and photometry tasks of the FITSH package (iPal 1201211 . 
As a hint, we used the MFC predictions for the target co¬ 
ordinates and then we performed individual centroid fit on 
each image and smoothed with polynomial regression for a 
better (much more precise and accurate) input for aperture 
photometry. Instrumental magnitudes were then extracted 
for both the reference stars and the target itself and then 
after applying the standard photometric transformations, 
we obtained the intrinsic g', r' and i' magnitudes. 

We searched for possible light curve variations using the 
most frequently sampled Sloan r’ band data (every second 
image was taken in Sloan r' while every fourth image was 
in g' and i'). In order to search for periodic variations in 
our data set, we fitted a function in a form of f{t) = a + b- 
sin(27rnt) -I- c • cos(27rnf) to the Sloan r' photometric series 
where t indicates the independent value (the time in this 
case). If the value of n is scanned in the appropriate domain 
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{n = 0.01... 15) with a proper stepsize (n = 0.01, that is 
^ 8 times smaller than the stepsize implied by the Nyquist 
criterion), then the parameters a, b and c can be obtained 
via a simple weighted linear least squares fit procedure. 
The unbiased values can then be compared with the 
reference value of Xq . This reference value is obtained when 
n is set to 0 and the error bars are scaled by a factor of 
1.51 to yield a Xo equivalent to the degrees of freedom. The 
difference between the Xo a-^d the x^ related to the adopted 
period tells the significance of the detection while various 
additional possible periods can also be checked and/or ruled 
out according to the difference between the respective x^ 
values. We found a significant variation (xo — = 24.2) 

with a corresponding amplitude of Ar' = 0.045 ±0.007 that 
has a frequency of n = 5.11±0.12d“^. The folded light 
curves are displayed in Fig. [21 The mean magnitudes of 
these observations were g) = 20.274 ± 0.013, r'^ = 19.519 ± 
0.009 and = 19.316 ± 0.013. 

We have to note here that due to the daily aliases, the 
peaks around n ± 1 d~^ are also remarkable and there is a 
non-negligible chance that one of these frequencies belong 
to the intrinsic rotation of the object. The peak at n = 6.11 
has a x^ value which is only larger than that of the main 
peak by 3.5. 

In general, minor bodies in the Solar System feature 
double-peaked light curves. Hence, the rotational frequency 
of 2013 AZgo is more likely rirot = n/2d“^, equivalent to a 
period of Prot = 9-39 ± 0.22 h. In order to test the signifi¬ 
cance of a double-peaked light curve solution, we coadded a 
sinusoidal component with half of the frequency to the pri¬ 
mary variations. The amplitude of this component is found 
to be 0.013±0.008mag. This is only a 1.7 -ct detection, how¬ 
ever, a good argument for confirming the assumption for an 
intrinsic rotation period of Ppot ~ 9.4 h. 

In addition, we repeated the photometric observations 
for 2013 AZeo in 2014 January 28 in g' and r' bands. 
The results of these photometric measurements yielded 
the Sloan magnitudes of g 2 = 19.71 ± 0.04 and = 
18.99 ± 0.03. During the first series of measurements (in 
2013 November), the geocentric and heliocentric distance 
of 2013 AZgo were Ai = 8.176 AU and ri = 8.244 AU, re¬ 
spectively, while in 2014 January 28, these distances were 
A 2 = 7.148 AU and ^2 = 8.114 AU. Based on these dis¬ 
tances, the expected change in the apparent brightness 
was 5[logiQ(r2A2) — logiQ(riAi)] = —0.326, however, the 
actual brightness changes were Ag' = —0.56 ± 0.04 and 
Ar' = —0.53 ± 0.03, whose mean is Am = —0.54 ± 0.03. 
Since the phase angle of 2013 AZeo was ai = 6.9° in 2013 
November 5 and a 2 = 1.5° in 2014 January 28, these 
values imply a phase correction factor of /3 = [(0.54 ± 
0.03) - 0.326]/(6.9 - 1.5) = 0.039 ± 0.006mag/deg. This 
is is rather good agreement with MFC observations. Based 
on the MFC observation database, the best-fit phase correc¬ 
tion parameter can also be derived, however, with a larger 
uncertainty: /3mpc = 0.040 ± 0.025 mag/deg. 

These parameters allowed us to derive the absolute 
brightness of the object 2013 AZgo in a manner described 
in the following. First, we employed simple Monte-Carlo 
run whose input were the observed Sloan brightnesses, the 
derived phase correction factor as well as the parameters 
and the respective uncertainties of the corresponding Sloan- 
UBVRl transformation equation (for conve rting g' and r' 
brightnesses to V, see [Jester et al.L l2005ll . This Monte- 
Carlo run yielded a value of Hy = 10.42 ± 0.07. Next, 


we checked the available photometric data series presented 
in the MFC database which yielded slightly fainter val¬ 
ues, namely Hv,mpc = 10.60 ± 0.15. In order to reflect 
MFC photometry in our derived absolute brightness value, 
we adopted the weighted mean value of these two values, 
namely Hy = 10.45 with a conservative uncertainty of 
±0.10 in the subsequent modelling. 

2.3. Reflectance spectrum 

In order to accurately compare the surface colour charac¬ 
terist ics of 2013 AZgo with other TNOs Isee lLacerda et all 
[ 20 TI . we obtained a low resolution spectrum using the 
Optical System for Imaging and Low Resolutio n Integrated 
Spectroscopy (OS IRIS) camera spectrograph (ICepa et al.L 
I 2 OOOI : ICepalDOl^ at the 10.4m Gran Telescopic Canarias 
(GTC), located at the El Roque de los Muchachos 
Observatory (ORM) in La Falma, Canary Islands, Spain. 
The OSIRIS instrument consists of a mosaic of two Marconi 
CCD detectors, each with 2048 x 4096 pixels and a total 
unvignetted field of view of 7.8' x 7.8', giving a plate scale 
of 0.127"/pixel. However, to increase the signal to noise for 
our observations we selected the 2x2 binning mode with a 
readout speed of 200 kHz (that has a gain of 0.95e“/ADU 
and a readout noise of 4.5e“), as corresponds with the 
standard operation mode of the instrument. A 300 s expo¬ 
sure time spectrum was obtained on January 28.17 (UTC), 
2014 at an airmass of A = 1.14 using the OSIRIS R300R 
grism that produces a dispersion of 7.74 A/pixel, covering 
the 4800 — 9000 A spectral range. A 1.5" slit width was used 
oriented at the parallactic angle. 

Spectroscopic reduction has been done using the stan¬ 
dard IRAF tasks. Images were initially bias and flat- 
held corrected, using lamp hats from the GTC Instrument 
Calibration Module. The two-dimensional spectra were 
then wavelength calibrated using Xe±Ne±HgAr lamps. 
After the wavelength calibration, sky background was sub¬ 
tracted and a one dimensional spectrum was extracted. To 
correct for telluric absorption and to obtain t he re lative 
rehectance, G2V star Landl02_1081 l|Landoltl . Il992f) was 
observed using the same spectral conhguration and at a 
similar airmass immediately after the Gentaur observation. 
The spectrum of the 2013 AZgo was then divided by that of 
Landl02_1081, and then normalized to unity at 0.55/rm. 

The derived spectrum is displayed in Fig. [3] Based on 
this spectrum, the slope parameter of this object is found 
to be S' = 13.4 ± 3.0 %(1000A)“^ by a linear ht across the 
interval 5000 — 9000 A. 

The measured photometric colours (g' — r' = 0.755 ± 
0.018 and 0.72 ±0.05 on 2013.11.04 and 2014.01.28, respec¬ 
tively) are in complete accordance with the derived spec¬ 
tral slope. The spectrum was normalized at 5500 A, just 
between the g' and r' bands . We can therefore write S' to 
equation (2) of jJewittl l|2002fl . if we write SDSS colours in¬ 
stead of Bessel ones, and set AA = 1480 A. This results in 
a synthetic colour index from spectral slope {g — r)synth = 
0.76 ± 0.04, in a perfect agreement with our photometry 
within the errors. 

2.4. Near-infrared photometry 

CCD observations of 2013 AZgo were obtained on 24 
September 2013 with the 4.2-m William Herschel Telescope 
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Fig. 4. Spectral energy distribution of 2013 AZeo in the far-infrared region, based on Herschel/PACS measurements. Left 
panel: far-infrared measurements superimposed are the best-fit NEATM curves with their respective uncertainty. Middle 
panel: TPM model curve for thermal inertia of 100 Jm“^K“^s“^/^, rotation period of 9.39 h and equator-on geometry. 
Right panel: the value of as the function of thermal inertia for pole-on and equator-on geometries. See text for further 
details. 


Table 3. Orbital and optical data for 2013 AZeo at the time 
of the Herschel observations. 


Quantity 

Symbol 

Value 

Heliocentric distance 

r 

8.702 AU 

Distance from Herschel 

A 

8.560 AU 

Phase angle 

a 

6? 6 

Absolute visual magnitude 

Hy 

10.45 ±0.10 


at La Palma Observatory, equipped with the LIRIS in¬ 
strument. LIRIS is a near-lR imager/spectrograph, which 
uses a Ik X Ik HAWAII detector with a field of view of 
4.27' X 4.27'. The number of exposures taken in different 
filters are: 5 x 30 s in T and J, 15 x 20 s in H, 15 x 13 s 
in CH4, and 180 x 20s exposures in Kg. Local comparison 
stars were selected from the 2MASS catalo gue and magni¬ 
tude t ransformation were applied following Hodgkin et al.l 


The result of the photometry is T = 18.66 ± 0.08, J = 
18.34±0.05,i7 = 18.00±0.06 and Kg = 17.72±0.10 where 
Y refers to the UKIDDS system, while JHKg are 2MASS 
magnitudes. Thus, 2013 AZgo exhibits almost exactly so¬ 
lar colour indices, with a slightly redder slope than solar. 
Namely Y - J = 0.32, J - H = 0.34 ±0.07 and H - K g = 
0.28±0.11 while according to lCasagr ande et HI (120121) an d 
estimating solar Y — J according to iHodgkin et ah ( 200911 , 
the respective solar colours are {Y — J)q = 0.235 ± 0.018 
(J - H)q = 0.286 ± 0.018 and {H - = 0.076 ± 

0.018. We note here that LIRIS is equipped with Mauna 
Kea Observa tories (MKO) system of J, H and K filters. 
According to iHodgkin et al.l (l2009rFI . the expected system¬ 
atic differences between 2MASS and LIRIS/MKO colours 
of 2013 AZgo is in the range of —0.010... ± 0.015, which 
is definitely smaller than the photometric uncertainties. 
This observation indicates a flat and featureless spectrum 
of 2013 AZgo: the slope is equivalent in the infrared and in 
the optical, being quite similar to dormant cometary nuclei. 

During the observation the heliocentric and geocentric 
distance of 2013 AZeo were 8.87 and 8.32 AU, respectively, 
indicating an absolute mid-IR brightness of J = 9.00 ±0.06 
without correcting for the solar phase angle. 


See their equations (6), (7) and (8) 


3. Thermal emission modelling 

3.1. Near-Earth Asteroid Thermal Model 

The basic physical properties such as albedo and diameter 
can be obtained by combining optical brightness data with 
thermal emission. Assuming an absolute optical brightness 
for a certain object, the higher the thermal emission, the 
smaller the actual albedo and hence the larger the diam¬ 
eter. Due to the observing strategies constrained by the 
spatial attitude of the Herschel, objects close to the eclip¬ 
tic are likely observable by Herschel during quadratures. 
Nevertheless, in quadratures these minor objects exhibit a 
large phase angle, hp rce Standard Thermal Model (STM 
iLebofskv et al.i r fT986ll might not be as accurate as it is de¬ 
sired. 2013 AZeo had a phase angle of a = 6 . 6 ° at the time of 
our Herschel/PACS observations. In order to have accurate 
estimates for larger phase angles, we empl oyed th e Near - 
Earth Asteroid Thermal Model (NEATM iHarrisL Il998ll : 
this model integrates the thermal emission for arbitrary 
viewing angles. Throughout the modelling we use the he¬ 
liocentric and geocentric distances at the instance of the 
Herschel/PACS measurements, namely Cheiio = 8.702 AU 
and rgeo = 8.560 AU. 

The diameter and albedo can then be deriv ed in a sim¬ 
ilar manner like in ou r earlier works (see e.g. iKiss et all 
l2013l iPal et all l2012ll . As an input for the fit procedures, 
we used the previously obtained thermal fluxes (see Table [21) 
and the absolute brightness Hy = 10.45±0.10 (derived ear¬ 
lier, see above). 

The absolute physical parameters (diameter, albedo and 
beaming parameter) have been obtained in a Monte-Carlo 
fashion. In each step, a Gaussian value were drawn for the 
four input values (three thermal fluxes and the absolute 
brightness Hy) and the model parameters were adjusted 
via a nonlinear Levenberg-Marquardt fit. A sufficiently long 
series of such steps yields the best fit values as well as the 
respective uncertainties and correlations. This procedure 
were performed in two iterations. Eirst, we let the value 
for the beaming parameter rj to be floated. This run re¬ 
sulted relatively high correlations between the parameters 
and a highly long-tailed distribution for 77 : we found that 
the mode for rj was 0.8 while the median is 2.6 and the 
uncertainties yielded by the lower and upper quartiles are 
This skewed distribution is due to the fact that 
beaming parameters cannot really be constrained if thermal 
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fluxes are not known for shorter wavelengths (i.e. shorter 
than the peak of the spectral energy distribution). Hence, in 
the next run we used ry as an input (instead of an adjusted 
variable) while its value was drawn uniformly between 0.8 
and 2.6. This domain is also in accordance with the possible 
p hysical domain of the beaming parameter (see also Fig. 4 
of iLellouch et al.l . l2013ll . The results of this second run were 
d = 62.3 ± 5.3 km, pv = 0.029 ± 0.006 while the beaming 
parameter can be written as rj = 1.7±0.9. The resulting 
albedo refers to a remarkably dark surface. The fluxes along 
with the best-fit NEATM model curve are shown in the left 
panel of Fig. S] 


3.2. Thermophysical Model 

In addition to the derivation of the NEATM parameters, we 
conducted an analysis of thermal emiss ion based on the as¬ 
teroid ther mophysical model (TPM, see lMiiller fc Lasrerrosl 
I1998L[200^ . The observational constraints employed by this 
model was the thermal fluxes (see Table. HD, the absolute 
magnitude of Hy = 10.45±0.10 (see earlier), the rotational 
period of 9.39 h as well as the actual geometry at the time 
of Herschel observations (see the values for phase angle, 
heliocentric and geocentric distances above). 

Our procedures have shown that the best- 
fit model occurs at high thermal inertia values. 
Assuming an equator-on geometry, a value for re¬ 
duced ± 1 corresponds to F > 10 Jm“^K“^s“^/^ 
(see also the right panel of Fig. [4]), however, the 
gradually decreasing form of the function x^(r) im¬ 
plies a lower limit of F > 50 Jm ^s The cor¬ 
responding values at F = 50 Jm“^K“^s“^/^ for ge¬ 
ometric albedo and diameter are py = 0.028 and 
d = 64.9 km, respectively. These values are also com¬ 
patible within uncertainties with the ones derived from 
NEATM analysis (see above). The spectral energy distri¬ 
bution provided by these TPM values are shown in the 
middle panel of Eig. SI This value for the thermal iner¬ 
tia is close to the values of 30 — 300 Jm“^K“^s“^/^ re- 


ported for comets llulia n. Samarasinha fc Beltonl. 20001 


ICampins fc Fernanded . 120001 iDavidsson et al.L 120131 see 
e.g.) as well as the value of 1 0 — 50 Jm~^K~^s~^ /^ for 
67P/Churyumov-Gerasimenko (|Gulkis et al.l . 1201 5ll . We 
note here that models either with thermal inertia values 
smaller than 50 Jm“^K“^s“^/^ or having an assumption for 
pole-on view underestimate the observed flux at 160 pm. 

Our findings for large preferred values of the 
beaming parameter p as well as for the thermal in¬ 
ertia F (even F > 10 Jm“^K“^s“^/^) ca n be compared 
with t he statistical expectations of iLellouch et al.| 
dml). By considering the small heliocentric distance of 
this object, both of th ese values are expect ed to be smaller 
(see Figs. 6 and 11 in ILellouch et al.l . l2013ll . 


4. The dynamics of 2013 AZgo 

2013 AZgo moves on a highly eccentric orbit, with perihe¬ 
lion between the orbits of Jupiter and Saturn. The best 
fit solution for the epoch of March 15, 2015 is shown in 
Table 01 

In order to assess the dynamical history, and po¬ 
tential future behaviour of 2013 AZ^n, w e follow a well- 
established route (see e.g. iHorner et al.l . l2004al lbl. 120101 


Table 4. The best-fit orbital solution (semi-major axis a, 
perihelion distance q, eccentricity e, inclination i, longi¬ 
tude of ascending node B, argument of perihelion w, mean 
anomaly M and perihelion date Tperi), and associated un¬ 
certainties, for 2013 AZqo, taken from the Minor Planet 
Center (http://www.minorplanetcenter.net) on March 15, 
2015. Note that the uncertainties of the orbital elements 
involved throughout the planning of the observations sig¬ 
nificantly larger than these due to the shorter arcs available 
at that time. 


a (AU) 

829.7 

q(AU) 

7.908098 ± 0.000014 

e 

0.990468 ± 0.000010 

i (deg) 

16.535760 ± 0.000011 

Q (deg) 

349.21122 ± 0.00002 

w (deg) 

158.14327 ± 0.00021 

M (deg) 

0.00876 

Tperi 

2456988.0641 ± 0.0032 

Epoch (JD) 

2457200.5 


I 2 OI 2 I : iKiss et ahl . 120131) . and used the Hybrid integra- 
tor within the n-body dynamics package MERCURY 
(|Chambersl . Il999l) , to follow the evolution of a swarm of test 
particles centered on the best-fit orbit for the object, in or¬ 
der to get a statistical overview of the object’s behaviour. 
A total of 91,125 test particles were created, distributed 
uniformly across the region of orbital element phase space 
within ±3crof the best-fit perihelion distance, q, eccentric¬ 
ity, e, and inclination, i. In this manner, we created a grid of 
45 X 45 X 45 test particles distributed in even steps across 
the ± 3(7 error ranges about the nominal best fit orbit in 
each of the three orbital elements studied. Each of these 
test particles was then followed in our integrations, with its 
orbit evolving under the gravitational influence of the giant 
planets Jupiter, Saturn, Uranus and Neptune, for a period 
of four billion years. Test particles were considered to have 
been ejected from the Solar System (and were therefore re¬ 
moved from the integrations) if they reached a barycentric 
distance of 10,000 AU. Similarly, any test particles that col¬ 
lided with one of the giant planets, or with the Sun, were 
removed from the simulations. Each time a test particle 
was removed in either of these manners, the time at which 
the removal occurred was recorded, allowing us to track the 
number of surviving test particles as a function of time. The 
results of our simulations are shown below, in Figs [5] and 01 

It is immediately apparent that the population of clones 
of 2013 AZgo is highly dynamically unstable, with 63.9% 
of the particles (58191 of 91125) being removed from the 
simulations within the first million years of the integrations, 
as a result of either ejection or collision with one of the giant 
planets or the Sun. Half of the test particles are ejected 
within the first 682 kyr of the integrations, revealing that 
the orbit of 2013 AZgo is more than two orders of magnitude 
more unstable th an that of the similar object 2012 DR 30 
(|Kiss et al.l . 120131) . 

The orbit of 2013 AZgo proves to be highly dynami¬ 
cally unstable on timescales of just a few hundred thousand 
years. Fully half of the test particles in our simulations were 
removed from the simulations within just 682 kyr, and al¬ 
most two-thirds were removed within the first million years. 
This extreme level of instability is not, however, that sur- 
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Logjo (Time Elapsed (yr)) 



Logjp (Time Elapsed (yr)) 


Fig. 5. The decay of our population of 91,125 clones of 
2013 AZgo as a function of the time elapsed in our integra¬ 
tions. The plots on the right show the same data as those 
on the left, but are plotted on a log/log scale. The upper 
panels show the decay of the population over the first bil¬ 
lion years of the four billion year integrations, whilst the 
lower panels show the decay over the first million years. 
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prising - 2013 AZeo passes through the descending node of 
its orbii |3 at essentially the same time it passes through per¬ 
ihelion, maximizing the likelihood that it will be perturbed 
by either Jupiter or Saturn. This extreme level of instabil- 
ity is typical of objec ts moying on Centaur like orbits (e.g. 
iHorner et al.l . I2004al f9l and suggests that 2013 AZgo may 
only recently haye been captured to its current orbit. This 
argument is supported by the fact that, ayeraged oyer our 
entire population of 91,125 test particles, the mean lifetime 
of 2013 AZeo is just 1.56 Myr. 

Given that 2013 AZgo exhibits such extreme instability, 
and may well be relatively pristine object, it is interest¬ 
ing to consider whether it will have experienced significant 
solar heating, and cometary activity, over its past history. 
As a result of our large dynamical dataset on the evolu¬ 
tion of 2013 AZgo, it is possible to determine the fraction 
of the population of clones that may one day evolve onto 
Earth-crossing orbits, and the fraction of the population 
that approach the Sun to within a given heliocentric dis¬ 
tance at some point in their lifetime. Since dynamical evolu¬ 
tion under the influence of gravity alone is a time-reversible 
process, we can use these values to estimate the probability 
that 2013 AZgo has moved on orbits that bring it within 
those heliocentric distances at some point in the past, be¬ 
fore being ejected to its current orbit. Due to the extreme 
instability exhibited by 2013 AZgo, we found that a rela¬ 
tively small number of the total population of clones were 
captured to Earth-crossing orbits through their lifetimes. 

® as can be seen in the elegant Java visualization 
of the object’s orbit at http://ssd.jpl.nasa.gov/sbdb.cgi? 
sstr=2013%20AZ60;orb=l;cov=0;log=0;cad=0#orb 


Fig. 6. The mean (upper) and median (lower) lifetimes of 
2013 AZgo, as a function of the initial perihelion distance, 
q, and eccentricity, e, of the orbit tested. The location of 
the best-fit orbital solution for 2013AZgo, as detailed in 
Table SI is shown by the hollow square at the center of 
the figure, with the Tier uncertainties on the perihelion 
distance and eccentricity denoted by the solid black lines 
that radiate from that box. Each coloured square in the fig¬ 
ures shows the mean (or median) lifetime of the 45 individ¬ 
ual runs carried out at that particular a-e location. Each 
of those 45 runs tested a different orbital inclination for 
2013 AZeo, evenly distributed across the ±3 (t uncertainty 
range on the best-fit orbital solution. As was the case wit h 
the high-eccentricity Centaur 2012 DR 30 (|Kiss et al.l . l2013ll . 
the stability of the orbit of 2013 AZgo does not vary signif¬ 
icantly across the range of perihelion distance and eccen¬ 
tricities tested in this work - a reflection of the relatively 
high precision with which the object’s orbit is known. 


Indeed, just 3805 of the 91125 test particles we studied (just 
4.2% of the population) became Earth-crossing at any point 
in our integrations, and the total fraction of the object’s 
lifetime spent as an Earth-crossing object (averaged across 
all 91125 test particles) was 0.12%. Our results for a variety 
of other perihelion distances are displayed in Table [SJ to¬ 
gether with estimates of the mean amount of time for which 
clones of 2013 AZgo exhibited perihelion distances smaller 
than the specified value. 
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Number of clones 

Percentage of total 
integration time 

Earth-crossing 
{q < 1.0616 AU) 

3805 

0.118 

q < 2 AU 

6005 

0.291 

q < 4AU 

12272 

0.329 

q < 6 AU 

27150 

2.06 


Table 5. The number of the 91,125 clones of 2013AZ6O 
simulated in this work that evolved to orbits with perihe¬ 
lion distances smaller than 2, 4, and 6 AU, and the number 
that evolved onto Earth-crossing orbits (following Horner 
et al., 2003). For each of these values, we also give the frac¬ 
tion of the total integration time, across all 91,125 clones, 
for which clones have perihelion distances within these lim¬ 
its. We note that this is the fraction of the time for which 
the perihelion distance was less than the stated amount 
and not the fraction of time the clones spend within that 
heliocentric distance. Even when moving on an orbit with 
perihelion within that of the Earth, a given clone will spend 
the vast majority of its time beyond that distance, and only 
a tiny fraction within it. 



Fig. 7. Slope parameter vs. albedo relations for 111 TNOs, 
including 2013 AZgo and 2012 DR 30 . Data (e xcept for these 
two la tter objects) have been taken from iLacerda et al.l 
(|2014ll . The b lue and red dot s indic ate the two major groups 
identified bv ILacerda et al.l (|2014ll . black points represent 
ambiguous objects (due to their large respective uncertain¬ 
ties) while green and yellow dots show the large bodies 
and Haumea-type surfaces, respectively. The isolated pur¬ 
ple square shows the place of 2013 AZgo at the very left 
side of the diagram. The other purple square indicates 
2012 DR 30 , just in between of the dark neutral (blue) and 
bright red (red) object groups. 


5. Discussion 

Since the orbit of 2013 AZgo is highly eccentric, and takes 
the object out to approximately 1950 AU, it is clear that 
it spends the vast majority of its orbit at large heliocen¬ 
tric distance. It is quite plausible that 2013 AZgo is a rel¬ 


atively recent entrant to the inner Solar System. Hence, it 
is interesting to consider how much time, cumulative over 
its entire history since it was first emplaced on a planet 
crossing orbit, has it spent at a heliocentric distance of less 
than 1, 10 or 100 AU. Again, we can take advantage of the 
large dynamical dataset available to us from our integra¬ 
tions to get a feel for the amount of time the object will 
have spent within these distances. Clearly, this is only an 
estimate (and implicitly assumes that, prior to its injection 
to a planet-crossing orbit, the object was well beyond the 
100 AU boundary - i.e. that it was injected from the inner or 
outer Oort cloud, rather than the trans-Neptunian region). 
Given that implicit assumption, we find that, on average, 
clones of 2013 AZgo spend just 6.68 years within 1 AU of the 
Sun, 4620 years within 10 AU of the Sun, and 273,000 years 
within 100 AU of the Sun. The time spent within 10 and 
100 AU is strongly biased by a few particularly long-lived 
clones, especially those that are captured onto Centaur-like 
orbits. We note that more than two-thirds of the clones 
(64904 objects) spent less than a thousand years within 
100 AU of the Sun, and 37493 (41.1%) spent less than one 
hundred years within 100 AU. Taken into considerations, 
our dynamical results suggest that 2013 AZgo has only re¬ 
cently been captured to its current planet-crossing orbit, 
and that it is quite likely that it is a relatively pristine 
object. Indeed, it seems highly probable that the surface 
of 2013 AZgo has experienced only minimal outgassing and 
loss of volatiles since being captured to a planet crossing 
orbit, and so it represents a particularly interesting object 
to target with further observations as it pulls away from 
the Sun following its recent perihelion passage. 

Outer Solar System objects can also be characte rized 
in a way recently put forward by ILacerda et al.l (|2014f) , us¬ 
ing their visual range colours and albedos. In this frame. 
Centaurs and trans-Neptunian objects form typically two 
clusters, a dark-neut ral and a bright-red one (see Fig. 2 in 
ILacerda et al.l . [20T3i . In this scheme, 2013 AZgo is located 
at the dark (very low albedo) edge of the dark-neutral clus¬ 
ter, see Fig. [T] 2013 AZgo is even darker than the object 
20 02 GZ 3 sj, the object w ith the lowest albedo in the sample 
of iDuffard et al.l (f2014f) . Objects with characteristics simi¬ 
lar to our target belong rather to “dead comets” or Jupiter 
family comets which are the end states of Centa urs and 
Oort cloud comets (Fig. 4 in ILacerda et al.l . l2014fl : in this 
sense 2013 AZeo is more similar to objects in the inner Solar 
System than those in the trans-Neptunian population. We 
also checked the distribu tion of the slope param eters of var¬ 
ious Centaurs based on iFornasier et all (120091) . Although 
in that work, a correlation between the slope parameters 
and orbital eccentricity were suspected (the higher the ec¬ 
centricity, the redder the object), the large eccentricity of 
2013 AZgo do not fit in this model since it has definitely 
lower slope parameter than the mean of that sample of 
Centaurs. 

While the dynamical analysis indicate that 2013 AZgo 
has recently been pulled from the Oort cloud, in the case 
of this object there is a much higher likelihood that it has 
spent a considerable time in the inner Solar System then 
e.g. in the case of 2012 DR 30 , which might just be in a tran- 
sitional phase bet ween the two main albedo-colour clusters 
(jKiss et al.l . 120131) . 
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